Origin of the highest energy cosmic rays 
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Introducing a simple Galactic wind model patterned after the solar wind we show that back-tracing the orbits 
of the highest energy cosmic events suggests that they may all come from the Virgo cluster, and so probably from 
the active radio galaxy M87. This confirms a long standing expectation. Those powerful radio galaxies that have 
their relativistic jets stuck in the interstellar medium of the host galaxy, such as 3C147, will then enable us to 
derive limits on the production of any new kind of particle, expected in some extensions of the standard model 
in particle physics. New data from HIRES will be crucial in testing the model proposed here. 



1. Introduction 

The origin of the highest energy particles ob- 
served in the universe continues to present a ma- 
jor challenge to physicists. These particles reach 
energies as high as 3 10^° eV. If protons, their 
flux is expected to drop sharply at 5 10^^ eV 
due to the interaction with the microwave back- 
ground, commonly referred to as the Greisen- 
Zatsepin-Kuzmin- or GZK-cutoff after its discov- 
erers. However, the number of particles known 
to be beyond 10^° eV continues to increase, with 
now 14 published |^,^,|4|, and a further 10 ex- 
pected from new observations with HIRES and a 
reanalysis of the Yakutsk data. 

There are three basic difflculties: 

First, we need to find a site that can produce 
particles at such energies [p^P] . There is only one 
class that has been argued to require protons at 
such energies in the source, namely radio galaxies 
with powerful jets and/or hot spots, such as M87 
m or Cyg A. M87 has been under suspicion to be 
the primary source for ultra high energy cosmic 
rays for a long time, e.g., ^,0. Gamma ray bursts 
are another possible class. 

Secondly, there is the additional difficulty of 
getting these particles to us, that is, overcom- 



ing the losses in the bath of the cosmological 
microwave background. That implies that the 
source should not be very much further than 20 
Mpc. 

The third difficulty is the explanation of the 
nearly isotropic distribution of the arrival direc- 
tions of these events. 

There are several ways out of these difficulties. 

Here we present the consequences of introduc- 
ing a magnetic Galactic wind in analogy to the 
solar wind. The magnetic field of the wind bends 
the particle orbits without adding substantial 
travel time. 

2. A model for a magnetic galactic wind 

It has long been expected that our Galaxy has 
a wind, e.g . ,H akin to the solar wind Recent 
modelling g shows that such winds can be quite 
fast, and ubiquituous. 

It seems plausible that this wind is powered by 
the combined action of cosmic rays and magnetic 
fields and starts in the hot phase of the interstellar 
medium seen in X-rays by ROSAT pl[| . 

Parker has shown that in a spherical wind 
the azimuthal part of the magnetic field quickly 
becomes dominant with Bfj, ^ sin0/r in polar 
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coordinates. 

Cosmic ray driving is similar to radiation driv- 
ing of winds in massive stars, and so [p^ magnetic 
fields can lead to an increase of the momentum of 
the wind. 

The data on the sign of the azimuthal compo- 
nent show that in the part of the sky most rel- 
evant for calculating orbits of energetic charged 
particles in our Galaxy, the field points to the 
direction of galactic longitude about 90 degrees 
That means immediately that positively 
charged particles traced backwards have their ori- 
gin above us, at high positive galactic latitudes. 

We adopt the simplest possible model. We as- 
sume that the magnetic field in the galactic wind 
has a dominant azimuthal component, and ignore 
all other components. We assume that this az- 
imuthal component has the same sign everywhere 
pl| . Most measures of magnetic field underesti- 
mate its strength. Therefore we will consider for 
reference a model which has a field strength near 
the Sun of 7 microGauss; this is close to the best 
estimates for the total local field |2|. The second 
parameter, the distance to which this wind ex- 
tends, is more uncertain: Our Galaxy dominates 
its near environment well past our neighbor, M31, 
the Andromeda galaxy, and might well extend its 
sphere of influence to half way to M81. Therefore 
we will adopt as outer the halo wind radius half 
the distance to M81, 1.5 Mpc. 

3. Tracing the path backwards 

To follow the particle trajectories in the Galac- 
tic halo we trace protons backwards, e.g., |22 1^, 
|l^ from their arrival direction at Earth. We use 
the 14 published cosmic ray events above 10^° 
eV, the list from Watson (included in Q) and 
the new list from AGASA [|^. There is a big 
uncertainty with the energy estimate of the high- 
est energy Yakutsk |l| event, which we therefore 
exclude from the present analysis, and hence we 
arrive at a final tally of 13 events used. 

The interesting result of these model calcula- 
tions is that the directions of all tracks point 
North. All events are consistent with arising orig- 
inally from Virgo A. Since these particles are as- 
sumed to be accelerated out of cosmic gas, about 



1/10 of all particles may be Helium nuclei with 
the same energy per particle. If the two highest 
energy events are in fact He nuclei, all 13 events 
point within 20 degrees of Virgo A. 

If Virgo A is indeed the acceleration site of 
the highest energy cosmic ray events, they all re- 
quire additional systematic bending at a ten to 
twenty degree level. Such bending could be easily 
accomodated within the plausible magnetic field 
strength within the supergalactic sheet from here 
to Virgo 
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How critical are our assumptions for these re- 
sults? 

The assumption of the symmetry of the mag- 
netic field above and below the Galactic disk is 
important. 

The value of the magnetic field, here adopted 
as 7 microGauss, for the wind near the Sun, is 
a key parameter. If the magnetic field were con- 
siderably weaker, the focussing would be largely 
removed. 

The scale of the Galactic wind here 1.5 Mpc, 
is not a critical parameter, since the calculations 
show that most of the bending happens within 
the first few 100 kpc. 

4. Discussion and implications 

This particular model can be tested in several 
ways; some of these tests are direct tests with 
data, and others take the form of a consistency 
check: 

The most important test is obviously with more 
data, and the HIRES data are soon to be released 
and will allow straight tests to be made. For in- 
stance, the apparent focussing is fully maintained 
at 5 10^9 eV. 

The very concept of a magnetic wind, driven by 
cosmic rays, but with an initial magnetic field as 
strong as that in the disk, needs to be examined 
more closely. 

There is one important consequence: If the 
model proposed here could be confirmed, then it 
would constitute strong evidence that all power- 
ful radiogalaxies produce high energy cosmic rays, 
and that they do this at a good fraction of their 
total power output. This then implies that com- 
pact radio galaxies u^M do provide a good test 
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bed for particle interactions, since they have a 
large screen of interstellar gas around the radio 
hot spots and jets as seen in mm- wavelength radio 
data There the paradigm is materialized of 
having a gigantic accelerator, and a beam dump. 
These radio galaxies may be used for particle in- 
teraction experiments in the sky. If a significant 
correlation in arrival direction between ultra high 
energy cosmic rays and this specific class of radio 
quasars could be confirmed ^ , then properties of 
new particles could be constrained. 

In summary, we propose here that a very simple 
model for a Galactic wind rather analoguous to 
the Solar wind, may allow particle orbits at 10'^" 
eV to be bent sufficiently to allow "super-GZK" 
particles to get here from M87, and also explain 
the apparent isotropy in arrival directions. 
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